Quantum Hall effect in bilayer and multilayer graphenes with finite gate voltage 
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We discuss the quantum Hall effect of bilayer graphene with finite gate voltage where the Fermi 
energy exceeds the interlayer hopping energy. We calculated magnetic susceptibility, diagonal and 
off-diagonal conductivities in finite-magnetic-field formalism, and observed crossover of integer quan- 
tum Hall effect from two independent monolayer type system to strongly coupled bilayer systems 
by changing the ratio of interlayer hopping energy and the gate voltage. We also discuss the case 
of multilayer systems with Bernal stacking. 
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Introduction — Among numbers of unusual physical 
properties of graphene (atomically thin graphite)^* 2 - 
study of the anomalous integer quantum Hall effect 
(QHE) was crowned with dramatic success. The Hall 
conductivity (per valley and per spin) is quantized as 

2 

a xy = — t~(ti + 1/2), n = 0, dbl,±2, ••• wheres a xy = 
— j^-n for usual two-dimensional (2D) electron gas. The 
difference 1/2 is explained theoretically, based on the 2D 
massless Dirac fermions where a Landau level is located 
at zero energy.— diSiil Moreover, QHE in bilayer system 

2 -, 

has also been observed as a xy = —^-n, n = ±1, ±2, • • • »± 
The characteristic feature of this result is that the step of 
a xy at the strong-magnetic-ficld limit is twice larger than 
other steps (see Fig. [TJ . This behavior was successfully 
explained by McCann and Fal'ko based on an effective 
Hamiltonian in 2 x 2 matrix fornix 

On the other hand, one of the special situations of 
graphene system which can not be realized in other Dirac 
fermion systems such as organic conductor a-(BEDT- 
TTF^Ig^ is that the Fermi energy is tunable parame- 
ter by the gate voltage. Although the theory by McCann 
and Fal'ko well describes QHE in sufficiently small Fermi- 
energy regions, their Hamiltonian is no longer valid when 
the gate voltage becomes greater than the energy gap be- 
tween two bands, which is the same order of the inter- 
layer hopping energy. This situation is recently realized 
experimentally^ For such cases, we need to discuss the 
QHE based on the Hamiltonian of the bilayer graphene 
in 4 x 4 matrix form [eq. ([T]) below]. 

The bilayer Hamiltonian is also important to discuss 
multilayer systems. For the Bernal stacking structure (it 
is also called AB or staggered stacking, and about 80 % of 
natural graphite falls into this category), the Hamiltonian 
of an A^-layer system given by 2N x 2N matrix can be 
block diagonalized into effective bilayer systems and a 
monolayer if N is oddi 10 ' 11 ' 12 ' 13 This decomposition is 
related to the Fourier modes of the wave function along 
the stacking direction^ In this paper, we discuss QHE 
of bilayer system with finite gate voltage based on the 



four band Hamiltonian (fT]), and also that of multilayer 
systems. 




(a) monolayer 



FIG. 1: (Color online) Schematic illustration of the anomalous 
integer quantum Hall effect (per valley and per spin) of (a) 
monolayer and (b) bilayer graphenes. The twice larger step at 
1/B — in (b) is due to the two- fold degeneracy of the Landau 
level at zero energy. The insets show dispersion relations of 
both systems. 

Eigenvalues and eigenstates — The Hamiltonian of the 
bilayer graphene in 4 x 4 matrix form is given by 



H = 
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where ir± = tt x ± m y with tv = p + eA/c is the momen- 
tum operator in a magnetic field V x A = (0,0, B). v 
and t are the Fermi velocity and interlayer hopping en- 
ergy, respectively. We have ignored the trigonal wrapping 
effect stems from the next-nearest-interlayer hopping for 
simplicity. For the system with zero-magnetic field H.q, 
dispersion relation is E )Ji (k) = s 2 (\/i 2 + (2vhk) 2 +sit)/2, 
where the label /i = (si,S2) specifies the outer and 
the inner bands (si = ±1), and positive and negative 
(,S2 = ±1) energies, respectively, so that t corresponds 
to an energy gap between two bands. Since the com- 
mutation relation between the momentum operators in 



cq. (U) is [ir±, ir T ] = ^f2eBh/c, there are correspondences 
with the creation and annihilation operators of the har- 
monic oscillator: n± — > \/2ja^ tt=f — > V2ja for eB ^ 0, 
where / = \J chj\eB\. Assuming that the wave function 
is given by linear combination of the number states of 
the harmonic oscillator \n), we obtain eigenvalues and 
cigenstates of eq. (U) as follows: The eigenvalues are 



tcring r = H/2t with r being the collision time of quasi- 
particles is given byi£ 
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E" = ^A" 



M = (si,s 2 ), 
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K = S2\ 
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where r = A t. For n > with A ^ 0, the cigenstates 

\/2hv ' 

are given by 
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The zero-energy state (A = 0) is doubly degenerate. For 
these states, we specify the quantum numbers n = and 
fi as follows, 
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For large interlayer hopping r 
expansion, the inner band si = 



(4) 
2 ~^> 1, using the Taylor 
1 is approximated as 



E-** = ± 



V2h 2 v 2 

t I 2 



y/n{n+l). 



(5) 



This is consistent with the result obtained by McCann 
and Fal'kofi Since the obtained eigenstates are given by 
number states like the nonrelativistic free fermion sys- 
tem, the degeneracy of each Landau level is also discussed 
in the same way that this is given by multiplicity of cen- 
ter of the coordinate: V/2irl 2 with V being volume of the 
system. 

Susceptibility — Before discussing the transport prop- 
erties, we derive the magnetic susceptibility based on the 
the finite-magnetic-field formalism to check the consis- 
tency with other formalisms. In fact, the Hamiltonian (TTJ) 
was introduced two decades ago by Safran motivated by 
graphite intercalation compounds^ He calculated mag- 
netic susceptibility based on the weak-magnetic field for- 
malism. An cxtcntion of this result with impurity seat- 



where f{x) = (c^ - ^' + l) -1 is the Fermi distribution 
function with /3 = l/k^T being inverse temperature. As 
discussed by Safran, in Y — > limit, the first term of 
the r.h.s of eq. (|6b[) gives diamagnetic logarithmic di- 
vergence near the zero Fermi energy |/x| — > 0, while the 
second term gives paramagnetic behavior around |/x| ~ t. 
Now let us calculate the magnetic susceptibility by the 
Landau quantization formalism. According to the func- 
tional integral method, the thermodynamic potential is 
given by 



Q(B) 
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Here w„ is Matsubara frequency of fermion including the 
chemical potential /x and effect of impurity scattering T 
as iuj„ = iuj„ + [fj, + i sga(u! n )T]/H. Applying the Euler- 
Maclaurin formula 
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g ( x )dx-±[g'(b)-g'(a)}, (8) 



to the second derivative of eq. ([7| : X — — h d 



V dW 



the 



s=o 



first (second) term of the r.h.s of cq. ([5]) gives the first 
(second) term of the r.h.s of eq. (|6b[) . Thus the result of 
this calculation coincides with that of the weak-magnetic 
field formalism. 

Conductivity — The conductivity is given by the Kubo 
formula as, 



m = 



Im Ilij (O 4- w?) 

hh : 



(9) 



where ELy(fi) 



ILij(n) -lhj(0) with {i,j} g {x,y}. 
The polarization function ILj- (Q) is given by the current- 
current correlation function, and obtained as the analyt- 
ical continuation of the Matsubara form: 
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where the matrix 7, is defined by 7 = Vfc7io/fi-, and 
Ei = Eu/h. \v m is Matsubara frequency of boson. 
We have ignored vertex corrections. After calculating 



the Matsubara frequency sums, and matrix elements 
((' ■'))((■■■ )) based on their symmetric (antisymmetric) 
properties for a xx (a xy ), the general expression of the 
conductivity (per valley and per spin) is obtained in the 
unified form as 



sgn(eB)a xy (ft) + ia xx (Q) 
e 2 2v 2 
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Here the matrix elements between Landau levels k and 
k + 1 remain for fc > 1, but other elements involving 
fe = states are complicated due to the degeneracy of 
zero energy state. We have defined 



(fi -> -n) (n^n) 



(12) 



and evaluated it analytically. 

Quantum Hall effect — For clean system T = and dc 
limit n - 0, X*" = [/(££) - f(Er)]/(E£ - E^ 2 . Then 
after straightforward calculations of cq. (fTTj) , expression 
of the Hall conductivity becomes a simple form, 



T X y = -sgn(eB) 
1 
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(13) 



Here, we have defined f(x) = f(x) — 1/2. In this deriva- 
tion it is important to keep the formula to satisfy the 
particle-hole symmetryi^i Moreover, assuming zero tem- 
perature and positive Fermi energy fj. > 0, we have 
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where [x]g means the integer part of x, and we have 
defined B = hv 2 eB/c. The first and the second terms 
of eq. (fT5|) correspond to contribution from the inner 
(,si = —1) and that from the outer (s± = +1) bands, 
respectively. As expected, the Hall conductivity is quan- 
tized as a unit of e 2 /h, reflecting that it is a topological 




FIG. 2: (Color online) Quantized Hall conductivity for T = 
[eq. (US}] versus inverse magnetic field 1/B for (a) bilayer sys- 
tems and (b) iV-layer systems (N — 1-5) with Bernal stacking 
and t/fj, = 0.5. 



number. Although, for the real QHE, the Anderson lo- 
calization due to the impurity scattering is essential, the 
present calculation with T = is useful to estimate pos- 
sible step structures of a xy . 

Results — The Hall conductivity a xy with T = versus 
inverse magnetic field 1/B at zero temperature [cq. (fT5"j) ] 
is shown in Fig. [Ha). Although the tunable parameter 
by gate voltage is fj, rather than t in experiment, we have 
shown the data as fixed t for later convenience. The 
results show clear crossover behavior from two indepen- 
dent monolayer like regime to strongly coupled bilayer 
regime as the ratio of t//j, is changed. Similar structures 
of a xy with inhomogeneous steps is also suggested to ap- 
pear in monolayer system with interaction between dif- 
ferent K points*^ In these data, we find that the slope of 
the Hall conductivity with respect to fj, 2 / B changes for 
\t/n\ > 1. This phenomenon is explained by rewriting 
inside of [• • • ]g in c q- (H3) as 



2\B\ 
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(16) 



In r.h.s of cq. HU), (■ • • ) is 1 + 0(\t/fi\) for \t/fi\ < 1, 
which means that the slope of conductivity is propor- 
tional to n 2 /B. On the other hand, for t//i > 1, the 
energy band related to the QHE is only the inner band 
(,si = —1), and (• • • ) becomes larger than 1 which means 
that the Landau levels strongly depend on the interlaycr 
hopping t. 

Next, in Fig.[3j we show 1/B dependence of the diago- 
nal a xx and the off-diagonal a xy conductivities with im- 
purity scattering T = 0.01/z. We have calculated dc limit 
of cq. (jlip analytically, except for summation over the 
Landau levels. In principle, this summation can also be 
performed analytically using the Poisson formula^ We 
find dip structure in the steps in the region \a xy \ > 5e 2 //i. 
The reason for this phenomenon is considered that the 
expression of a xy implies that of a xx , and contribution 
of <r xx becomes large as the magnetic field is decreased. 
Actually, in the classical result for nonrelativistic free 
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FIG. 3: (Color online) Crossover of quantum Hall effect from two independent monolayer like regime to strongly coupled bilayer 
regime, observed as changing the ratio of interlayer hopping energy t and the gate voltage (Fermi energy /i). The red and green 
lines denote off-diagonal (<J xy ) and diagonal ia xx ) conductivities with T/fi — 0.01, respectively. The blue lines indicate a xy 

with r = 0. 



fermions, the expression is a 



xy 



■ -^—a Tr where : 



and u> c oc eB are the electron density and the cyclotron 
frequency, respectively. This relation is also obtained in 
the analytical expression of the quantum case of 2D free 
fermions with constant TJ£ 

Finally, we consider multilayer systems. Since Hamil- 
tonians of multilayer systems with Bernal stacking are 
known to be decomposed into [N/2]g bilayer systems 
with effective hopping t* = t sin j^rrn I TO = — (-^ ~ 
1),—(N — 3), • • • , TV — 1], and a monolayer if N is 



odd j 10 i 1:L ' 12 i 13 the Hall conductivity can be calculated 
combining the above results. For example, the Hall con- 
ductivity with r = of iV-layer system (N = 1-5) ob- 
tained by eq. (fl"5|) and this matrix decomposition tech- 
nique with t/fj, = 0.5 is shown in Fig. E^b). Although 
QHE of multilayer systems with general closed packed 
stacking structure where Bernal and rhombohcdral stack- 
ing are mixed is discussed by Min and MacDonald using 
the matrix decomposition technique and the 2x2 effective 
Hamiltonians^ their results are limited to the small gate 
voltage regions where only the inner band contributes to 



Summary — We have discussed physical properties 
of bilayer graphene in finite magnetic field and finite 
gate voltage based on the Hamiltonian with four energy 
bands. We have checked the consistency between the 
weak- and the finite-magnetic field formalisms calculat- 
ing the magnetic susceptibility. Then general formula of 
conductivity in this systems has been obtained. We have 
observed crossover of integer quantum Hall effect from 
two independent monolayer system to strongly coupled 
bilayer systems by changing the ratio of interlayer hop- 
ping energy and the gate voltage. 
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After this paper has almost been completed, we be- 
came aware of the paper by M. Koshino and T. Ando, 
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4x4 Hamiltonian as ours to discuss the optical properties 
at zero gate voltage. 
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